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Abstract. We present an active tile assembly model which extends Winfree's abstract tile assem- 
bly model [31] to tiles that are capable of transmitting and receiving binding site activation signals. 
In addition, we introduce a mathematical framework to address self-similarity and recursion within 
the model. The model is applied to show a recursive assembly of an archetypal self-similar aperiodic 
tiling known as the L-shape tiling. 



1. Introduction 

DNA self-assembly is based on an inherent property of its chemical structure: a given nucleic 
acid sequence binds to its perfect Watson-Crick complement. Partially complementary sequences 
may also bind but these bonds are weaker. Moreover, DNA can assume a variety of rigid forms 
other than the familiar linear double helix. Synthetic molecules have been designed and shown to 
assemble into branched species |1H [30] , and more complex species that entail the lateral fusion of 
DNA double helices ^26j, such as DNA double crossover (DX) molecules [6j, triple crossover (TX) 
molecules [13] and paranemic crossover (PX) molecules. DX and TX molecules have been used as 
tiles and building blocks for large nanoscale arrays |3H l32] . 

DNA-based tile arrays can carry out computation: problems can be encoded in single-stranded 
DNA portions at the edges of the tiles (the "sticky ends" ) and the computation can be carried out by 
their self-assembly. Winfree introduced the abstract tile assembly model (aTAM) [31j and showed 
that 2D self-assembled arrays made of DX or TX DNA tiles can simulate the dynamics of a bounded 
one dimensional cellular automaton and so are potentially capable of performing computations as a 
Universal Turing machine. Since then, several successful experiments have confirmed the possibility 
of computation by array-like DNA self-assembly: binary addition (simulation of XOR) using TX 
molecules as tiles [16|, Sierpinski triangle assembly [23^ 17]. and a binary counter [2j by DX molecules. 
Transducer simulations with programmed inputs by TX DNA molecules have also been reported 

Si- 

DNA-based tiles also provide a way to obtain new materials both directly, as the array is a 
physical structure in itself, and indirectly, by allowing the tiles to contain markers or scaffolding 
which could direct the assembly of other molecular structures on top. Two-dimensional periodic 
tilings form crystal lattices and have been studied the most to date (see |36]). It is the aperiodic 
tilings which can carry out computation. In this paper we show a model for assembling self-similar 
aperiodic motifs: these, if materialized, could yield quasi-crystallographic structures, which are 
quite rare in nature and in general quite difficult to design and assemble. Our model is an extension 
of the abstract tile assembly model to tiles that are capable of transmitting and receiving binding 
site activation signals. One of the well-known aperiodic set constructions consists of modified square 
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tiles due to Robinson [22]. A tile set of just four "active" tiles utilizing the idea of signaling has 
been proposed in |18] for a set of Robinson tiles which are capable of hierarchically assembling a 
self-similar square pattern. In this paper, we propose a theoretical model for active tile assembly 
and a formal definition of self-similar and recursive assembly. We apply our model to another 
well-known aperiodic self-similar tiling, the L-shape tiling (also known as the chair tiling [8]). 

But why augment the existing model with signaling and how? In what has become the standard 
tile assembly model, the "sticky ends of the tile are formed by single stranded DNA containing a 
sequence of unpaired bases which can bind to a strand of complementary bases (on another tile). 
One important problem in the assembly process is error control: how to prevent "tiles" which have 
partial base matches from binding in the "wrong" place in the lattice. At high enough temperatures, 
weaker bonds are not stable; but, unfortunately, they can last long enough for the tile to form new 
and stable bonds in an incorrect position. One possible way to control the errors is by staging the 
assembly: by only allowing certain tiles in solution at a certain point, then adding other tiles in 
stages after a sufficient amount of time has passed for most correct bonds to form. A theoretical 
analysis of staged assembly is given in Demaine, et al. [Ij. 

Signaling is another possible solution, proposed in different variants in [15j and |18j . The idea 
is to keep certain binding sites dormant until a "signal" is received to activate them, which allows 
indirect control over the assembly process. The existing DNA signaling mechanisms rely on strand 
displacement techniques (see [35]) and have already been applied in a variety of ways. For example, 
devices whose activity is controlled by DNA strands have been produced in [35l |33j : these utilized 
DNA "fuel" strands, which were also recently used to construct a large number of enzyme-free logic 
gates [25l [20l |2T| . Additionally, structures that can perform simple "walking" (based on strand 
displacement) on an arranged platform have been reported in [271 ISl EH]. The three most recent 
walking devices |10| [T71 [Ti] showed a significant robustness, allowing the directions and the actions 
of the walker to be guided by a sequence of strand displacements incorporated within the walking 
platform. 

These mechanisms allow one to construct DNA based tiles with sticky ends that begin in an 
inactive state - bound to a protective strand which can later be displaced by a signal initiated 
by some other tile, resulting in activation. A 
particular DNA structure that can be used as 
a tile, an origami cross shape [36], is shown in 
Figure [l] (DNA origami is a method of creat- 
ing arbitrary shapes from DNA by folding a 
single strand into the desired shape and us- 
ing shorter "helper" strands to bind it rigidly 
together [S]). The size of this tile is about 
100 X 100 nm, which allows enough space to ac- 
commodate signaling, unlike the DX and TX 
molecules which may and have been used for 
computation [31] but are not large enough to 
incorporate the existing signaling mechanisms. 

We present the formal active assembly model 
in the next section, using a generalization of Figure 1. A DNA origami tile [55] , 

Winfree's aTAM known as the two-handed as- 
sembly model (2HAM) [1] (a part of what the authors of call a "multiple tile model"). In it. 
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the tile types are unique squares defined as 4-tuples {an,cre,as,(7w) representing the "glues" at 
the north, east, south, and west edges, each glue having an associated strength in the form of a 
non-negative integer. Tile assemblies are produced when smaller tile assemblies bind together. The 
"two" in 2HAM refers to the number of tile assemblies which may bind at a time. A bond is formed 
between adjacent tiles of two assemblies if their adjacent edges have matching glue labels and has 
the corresponding strength. Every tile assembly is a mapping of the integer lattice into a set of tile 
types and the bonds between all tiles are represented by a binding graph (the vertices are the tiles 
and edges are the connections among them). 

The formal model proposed in this paper is that of an active programmable self-assembly which 
is able to build up arrays hierarchically. In such a system, the signaling capabilities of the tiles allow 
them to assume multiple identities throughout the assembly process. Moreover, a tile is allowed to 
have more than one type of "glue" on any given edge and these may be initially inactive (unable 
to bind) until they receive an activation signal. This approach allows for step by step assembly in 
a potentially controlled and robust way without actually staging the process. 

We note that a formal description of the same process of "active assembly" allowing tile glues 
to be deactivated and tile assemblies to break apart was introduced in [19]. Although also based 
on the 2HAM, it differs in notation from our model and includes an added notion of "pending 
actions." The latter allows signal transmission to be delayed by an arbitrary amount of time in 
order to better approximate the physical interactions among actual DNA tiles. Both deactivation 
and "pending actions" are compatible with the model presented here, and we discuss the specifics 
of their incorporation into our model in the concluding chapter. 

The work presented here suggests a potentially feasible model for recursion in molecular self- 
assembly; and, through the use of signaling techniques, it has the potential to build self-similar 
structures hierarchically. In Section 2, we formally define active tiles and describe the hierarchical 
assembly in terms of a tile modification function. We also formally define the notions of recursive 
assembly and self-similarity. In Section 3, we apply the model to an archetypal example of a self- 
similar aperiodic tiling - the L-shape tiling p] - and prove that the designed tiling system recursively 
assembles a non-periodic self-similar L-shape based array, providing a tiling of the plane. We end 
the paper with some concluding remarks about the presented model and future avenues of research. 

2. Model Description 

In this section, we give the formal description of our Active Tile Assembly Model. In the 
first subsection, we define active tiles while in the second, we present binding graphs and tile 
assemblies. The dynamics of signaling within the assembly process are represented through the tile 
modification function which we describe in the third subsection. In the last two subsections we use 
the tile modification function to construct hierarchical tile assembly sets and to define an active 
tile assembly system as well as the properties of recursion and self-similarity. 

2.1. Basic Definitions. Let S+ be a finite set and define the complementary set = {— c| c G 
We call S = 11+ U the set of all labels. We use |c| to denote c if c S 11+ and to denote — c 
if c E S~. We also have — (— c) = c. We assume that S is fixed throughout our discussion. 

Definition 2.1. The strength function is defined to be 

s : S+ ^ Z+ where Z+ = {x G Z| x > 0} , 

and for c G S, s(|c|) is called the strength of c. 
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The strength of a label refers to the strength of the bond between it and its complementary label. 
Definition 2.2. Define the set of signals: 

t = jc^l c e i,j G {+y,+x, -y,-x,0} , j / o| , 

where cP- for i ^ means that c can be transmitted from the direction i in the direction j, whereas 
Cq means that c can be initiated in the j direction and we refer to it as an initiation signal. 

The direction notation in this paper follows the axial directions of the xy-plane. Although it is 
a departure from the conventional notation of North, East, South, and West p], we believe that 
this notation gives a more natural correspondence between the parts of a tile and its coordinate 
position in the plane and simplifies the description of signal transmission between adjacent tiles 
(e.g., we can say that a signal will "connect" to the signal of the form c^^; here, the notation j 
and —j indicates a more direct relationship between the corresponding sides than, say, and S). 

Definition 2.3. A tile t is a 4-tuple of tile sides, t = {t^y,t^x,t-y,t-x), such that each tile side ti, 
i G {+y, +x, —y, — x}, is an ordered pair of sets of labels: 

ti = {A, I) G X 

where CP(S) denotes the power set of S. We refer to A as the set of active labels and to / as the 
set of inactive labels. 

+y 



-X 




a) -y b) 

Figure 2. a) A tile t — {t^y,t+x,t-y,t_x) placed on the coordinate plane. Solid lines 
represent active labels and dashed lines represent inactive labels. Subscripts on tile sides 
indicate directions normal to those tile sides, b) A tile with complementary labels marked 
by +/-. See Example [2!4| 

Example 2.4. Refer to Figure [2] for an example tile t. Let = {r,b,y, g}, corresponding to the 
red, blue, yellow, and green labels. Then S~ = {— r, —6, —y, —g}, so S = {r, b, y, g, — r, —b, —y, —g}, 
and 



(W,W) 



(0,{-r,y,6,-5}) 



({-5} 



t. 



(W,0) 



are the tile sides for t = {t^y, t^^^t-y, t 



Definition 2.5. An active tile is an ordered triple r = {t,A,§>) where i is a tile and C S are 
the sets of activation signals and transmission signals respectively, such that 

(1) iorta = {A,I) 

a) if c G A, then — c ^ A and c, — c ^ / 

b) if d G /, then -d ^ I. 

(2) Cq ^ A for any a G {+y, +x, —y, —x} 
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Condition (1) states that a tile side cannot contain redundant or complementary labels; condition 
(2) requires that the set of activation signals does not contain initiation signals. 

Condition (1) comes from practical considerations of actual DNA dynamics: if a tile side contains 
active complementary ends, they may bind to each other instead of binding to complementary 
strands on another tile. Since we suppose that a signal for a given label would also activate the 
corresponding complementary label, we forbid complementary labels from co-existing even in the 
inactive state. Condition (2) eliminates superfluous signaling (i.e. a tile which activates one of its 
own labels is equivalent to a tile with that label already active) . See Fig jsja for an example of a tile 
which violates these conditions. 
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Initiation Signal 



Activation Signal 




Figure 3. a) An active tile obtained by adding signaling to tile t. See Example 



2.6 



-y^-g} 



b) 

Example of failed conditions (1), (2a) and (2b) in the active tile definition (Def. 2.5). This 
figure establishes a convention for depicting active tiles in the rest of the paper. 

Example 2.6. Refer to Figures [2] and [3^) for the sample tile with S = {r,b,y,g, —r, —b, 
and t = {t+y,t+x,t-y,t-x) as in Example |2.4[ The signaling shown in Figure [s] is: 

A = {rXl, h% yl:} 

and the corresponding active tile is r = (t,yi,S). ■ 

2.2. Tile Assemblies and the Associated Binding Graphs. Let T be any set of active tiles 
over an alphabet S with a strength function s and consider a partial mapping: 

72 



Q : 



1. 



We call such a mapping a configuration and define its associated binding graph to be the weighted 
graph Ga 
such that 



graph Ga = {V, E), with V C Z'^ such that (i, j) G ^ if and only if a{{i,j)) is defined, and E C (^) 



E — EfiU E^ 
Eh = {{{i,j),ii + l,j)}\ {i,j),{i + l,j)GV} 
Ev = {{{hi) , {hj + 1)} I ihj), + 1) G n • 
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The weight function for G is defined in the fohowing way. Let e = {v,w} = + 1)} be 

a vertical edge, and let 

a{v) = = (^^yl^S''), a{w) =T^ = {t^,X",$^) 

with 

J.II f^V J.V J.V \ J.W fj.W J.W J.W J.W \ 

~ \^+yT ''-y ''-xJi ~ V''+j/5 '"-x/' 

(where t^y = (A^yjI^y) and similarly for the rest). We consider the intersection 

where —A^y = {— c| c S ^-y}- Then we define the weight of the (vertical) edge e to be 

\ otherwise. 

The definition for the weight of horizontal edges is analogous. 

In other words, the vertices of the binding graph are formed by points in the integer lattice which 
we presume are occupied by tiles. There are horizontal and vertical edges indicating horizontal and 
vertical adjacencies. Each edge is assigned a weight corresponding to the sum of the strengths 
of the complementary active labels of the adjacent tile sides. We say that adjacent tiles whose 
adjacent sides contain complementary labels c and — c contain a bond of strength s(|c|). 

With this definition of a binding graph in hand, we proceed to make the following definitions. 

Definition 2.7. Given a set of active tiles T and an integer > called the temperature parameter, 
a tile assembly instance for is a configuration a where either 

(a) a is defined for only a single point in 1? or 

(b) its associated graph Ga = iV, E) satisfies the following two conditions: 

(i) G is connected, and 

(ii) the sum of the weights of the edges in any edge culj^of G is greater than or equal to 9. 

Observe that condition (b) in the definition, part (ii) of which is generally referred to as 9- 
stability, forces the tile assembly instance to be composed of at least two tiles; we refer to a in the 
special case of (a) as a unit tile instance. In literature, the temperature parameter 9 is usually set 
equal to 2. 

We use the terms tile assembly and unit tile to denote the equivalence classes [a] of the respective 
instances under translatior(3 

a ~ /? if and only if 3k, I £ Z such that a{i,j) = f3{i + k,j + I) G Z. 

We adopt the convention that the class representative (with respect to translation) of a tile 
assembly [a] is the tile assembly instance a' G [a] such that 

min |i | a'{i,k) is defined for some 

and 

min |j I a'{k,j) is defined for some A; G Z| = 0; 
that is, the tile assembly instance is positioned with the left-most tile on the y-axis and the bottom- 
most tile on the x-axis. 

-'^An edge cut of a graph G — {V, E) is a set of all edges with one vertex in U and the other in V\U for some 
[/ C V. 

o 

We do not allow rotation explicitly in the model. To account for rotation, we later require that each tile set 
contain "rotated" copies of every tile. 
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2.3. The Tile Modification Function. In this section we define the tile modification function / 

which simulates signal transmission and binding site activation for each tile assembly. This function 
allows each tile within the assembly to change its identity according to the signals it "sends to" 
and "receives from" its neighbors. 

Let a be a tile assembly instance and denote a{z), the tile associated with z, by az- For z 
and az = T = {t, A, S) G 7, let 

'^+2/ = «^+(0,l) 

T+x = a^+(i,o) 

T-y = a^+(0,-l) 

r-x = a^+(-i,o) 

define the neighbors of r on the respective sides. If Tj is not defined for some i G {+y, +x, —y, —x} 
write Ti = e, otherwise write r+y = ^1+^,8+^) (similarly for the other three neighbors). 
Let N = {(0, 1), (0, —1), (1, 0), (—1, 0)}. A function / : — >■ T is a tile modification function if 

/{cHz+n) = f{T , T+y, T+x, T-y, T-x) 

= a'z = {t',A',S') 

transforms the active tile r = into another active tile a'^ with the following three rules. 
Let t = {t+y,t+x,t-y,t-x) and t' = (t+j^, t+j., 

(1) Tile Side Modification 

If ti = {A,I), i e {+y, +x, -y, -x}, then 

t[ = {AuC, I\{CUD)), 

where 

C = |c G /| 3|c|*- G A and |c|q-' G S-' for some j G {+y, +x, —y, — x}| 

D = |c G /| $\c\^j G A for any j G {+y, +x, —y, — .t}} 

Informally, if the tile side contains an inactive label c with a corresponding activation signal 
and an adjacent tile contains a corresponding initiation signal, then the label c becomes 
active. Inactive labels which can never be activated are removed from the tile. 

(2) Activation Signal Modification 
For i, j G {+y, +x, -y, -x}, 

A A\ ^removed 

■^removed = |cj E A\ Ti ^ E and either 3cq * G or $c'^^ G for any k G {+y, +x, —y, — a;}| . 

Informally, when the tile receives a signal, its corresponding activation signal is used and 

cannot be used again, so it is removed. Activation signals that can never be activated due 
to being adjacent to a tile which has no possibility of sending the corresponding signal are 
also removed. 

(3) Transmission Signal Modification 
For i, j G {+y, +x, -y, -x}, 

S' = (§ U Sadded) \ ^removed 

Sadded = |c^| 3fc G {+?/, +x, -y, -x} such that G S and Cq G S'^j 
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"^removed 

5i U ^2 U S3 

Si = G §1 Tj / e and either 3cq * G 5* or ^c^* G S* for any k £ {+y, +x, —y, — x}| 
'S'2 = |4 ^ S| '^i 7^ ^ ^ctj G S-' U yi-' for any k G {+y, +x, -y, 

53 = {c^o e S| r,- / . 

Informally, when the tile receives a signal, its transmission signal is replaced with an ini- 
tiation signal. Transmission signals which do not have matching signals in adjacent tiles, 
i.e., transmission signals that cannot be received or transmitted, are removed since they 
cannot perform any function. Once a tile is adjacent to another tile, all initiation signals 
in that direction are removed, since they are either immediately transmitted or can never 
be transmitted. 

We can extend / to a tile assembly instance by applying / simultaneously to every active tile in 
the tile assembly instance: 

/(a) = a' if and only if a'^ = f{az+N) Vz G Z^. 

We write f{[a]) = [a'] = [/(a)] to denote the image of the tile assembly under /. 

Example 2.8. An example of / acting on an assembly of four tiles is shown in Figure |4j We will 
focus on the tile T = {{t^y,t-\-x,t-y, t-x),A, S) in the center of the assembly. Let r, b, g, y correspond 
to the red, blue, green, and yellow labels, respectively, in the figure. For side modification, only 
the labels in t+x = {^li^id}) are affected. In the first iteration, we have C = {g} and D = 0, 
because gty G /I is an activation signal for the green label and the tile T-y contains an initiation 
signal g^^ , so we end up with t'_^_^ = {{g},{r}). For the activation signals, we begin with A = 
{r^yyg-y} = -Aremoved becausc g^y was "activated" by gQ^ in the bottom tile and r^y does not 
have a corresponding transmission signal in the top tile T+y So, A' = 0. Finally, in the transmission 
signal modification, § = {63 y^J"^, ^, 5^^, 6^^}. Sadded = {^o^} since is adjacent to b^^ in 
T^y For the same reason b_^.y G Si. Si does not contain g^y because it is adjacent to g_^y in the 
T^y', however, g^y G S'2 because there is no corresponding transmission or activation signal in T_y; 
and, S3 = {bQ^,yQ^} because T has neighbors on the —x,—y sides but not on the +x side, so 
§' = {y+^b,y}. 

In the second iteration, C = and D = {r} because there is no longer an activation signal for r 
in A', so t'^^ = {{g},^)- A' = ^ = A" since activation signals can only be removed, never added, 
and S" = {y^^} because T has a neighbor on the —y side and therefore 6q ^ is removed. ■ 
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Figure 4. The tile modification function acting on a tile assembly. Two iterations of the 
tile modification function complete the modification of the assembly. Refer to Example 2.8 
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2.4. Hierarchical Tile Assembly Sets and the Active Tile Assembly System. We use the 

term "complete tile assembly" to refer to a tile assembly which has reached a quiescent state under 
iterated application of the tile modification function. 

Definition 2.9. A complete tile assembly instance for a is a tile assembly instance denoted a, such 
that 

a = /"(a) for some n and f{a) = a. 

As noted, [a] is the equivalence class of a under translation, so we define the complete tile assembly 
obtained from [a] as [a]. 

In figure |4j the tile assembly on the far right is complete. 

A seed set of unit tile representatives, denoted Tq, is a set of unit tiles that is closed under 90° 
rotations: for every a G Tq, if a' (0,0) can be obtained from a(0,0) through a cyclic permutation 
of the direction coordinates +x, —y, —x) in the tile description, then a' £ Tq. 

Fix a temperature 6. We define Tj, i > 1, recursively as follows: a £ Tj, if and only if a G 'Ji-i 
or a = a' where a' is a tile assembly representative composed of two tile assembly instances, f3i 
and f32- 

'( ) - I z £ dom(/3i) 

I f32iz) ifzGdom(/32) 

such that the domains of /3i and (32 are disjoint and such that there exist f3[,f32 G Tj-i with 
= [l3[] and [/32] = [/32]- In other words, a is in Tj if it is the completed instance of a composition 
of two tile assembly instances whose class representatives are in Tj_i. 
We call 7n an active supertile set at stage n. Observe that li-i C Tj. 

Definition 2.10. An Active Tile Assembly System (ATAS) is an ordered triple of a seed set, a 
strength function, and a temperature parameter: (To,s,^), 9 G Z^. 

2.5. Recursive Assembly. Now we proceed to define the recursive assembly of tile types and 
provide a bridge between the tiles of tiling theory and the unit squares of the TAM. We capitalize 
"Tiles" in what follows to emphasize that these are not necessarily unit squares. Intuitively, a 
Tile type could be thought of as a shape that consists of a fixed number of regions, each of which 
consists of subregions that all "look" the same (i.e., a region is a collection of instances of one kind 
of tile assembly) . The recursive assembly allows the number of subregions in each region to increase 
without changing the region's "kind," allowing the Tiles to grow into larger and larger versions of 
themselves. 

Continuing formally, let Too denote the set of all tile assembly representatives producible by an 

oo 

ATAS, i.e.. Too = |^ T„,- Then a Tile type is a function 

n=0 

r : Z+ ^ Too 

and T{i) is a Tile type T at level i. Furthermore, a Tile type T is self-similar if dom(T{i + 1)) is 
geometrically similar to dom(T(^)) for all ^ > 0. 
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A finite set = {Ti, . . . , T„} of Tile types is said to be recursive if for each ^ > 0, Tj(£) satisfies 
the following condition^ 

i=i 

where p{i) is the number of regions in the ilh Tile type, and Rij{i) is a configuration (but not 
necessarily a tile assembly instance), referred to as the jth region in the ith Tile type (at level £. 
We note that for each Tile type i, the number of regions p{i) does not depend on the level i. 

We require that the regions are unions of tile assembly instances whose number increases with 
the level. The regions Rij therefore satisfy 

dom (i?,,,, {£)) n dom {Ri^j, (^)) = for ji / 

and 

= U 

k=l 

where ijij : N — t- N is a non-decreasing function with rjij{i) specifying the number of subregions 
r^j{i) in the jth region of the ith Tile type at level £. 

The subregions r^-{£), rfj{i) must be assembly instances with non-overlapping domains for k ^ k' 
and further, VA;, k' > 0, 

(1) ye [r^^ii)] = \rf;^{i)] 

(2) (a) W,m [r^jie)] = [rfj{m)] or 

(b) Bio, q such that W K^^.(£)] = [Tg{eo + i)]. 
That is, all of the subregions in a given region must be instances of the same tile assembly for a 
fixed £ (condition 1) and also either be instances of the same tile assembly across all i or be Tile 
types (conditions 2 (a) and (b) respectively). Note that q and £o are specific to the jth region of 
the ith Tile type. 

Remark. Observe that condition 2(b) allows a Tile type i at level £ to contain a Tile type j ^ i at 
level i' as a subregion even if i' > i, since Ti{£) need not appear before Tj{i') does. One or more 
levels of some Tile types may assemble before the assembly of the first level of other Tile types 
begins. 

Thus, a Tile type in a recursive set is a sequence of tile assembly representatives indexed by level, 
each of which is the union of a number of configurations (regions) that is fixed for the Tile type and 
independent of level. Each of these regions, in turn, is the union of tile assembly instances that all 
belong to the same class but whose number is allowed to increase as the level increases. The class 
of instances corresponding to a given region, on the other hand, is either fixed and independent of 
the level or else encompasses a single Tile type. If the region R corresponds to a single Tile type 
T, then R at level 0, R{0), consists of translated copies of T{i') for some £' > and, in general, 
R{i) consists of translated copies of T{i' + i). 

Definition 2.11. A subset ©' = {Tj^, . . . ,Ti^} of a recursive set of Tile types © = {Ti, . . . ,T„} 
where each Tj^, is self-similar is called strongly self-similar if for every region Ri^j, k = 1, . . . ,m, 
there exists a Tile type Tj^, G O' such that for all I, dom(iij^j' (£)) is geometrically similar to 
doui{Ti^,). 

■^By union we mean the piecewise joining of the indicated configurations (recah that a configuration is a mapping 
from to the set of possible active tiles). 
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In other words, a subset G' of a recursive set Q is strongly self-similar if every region of each 
Tile type in 0' is geometrically similar to some Tile type in 0' (specific to the region). 

Remark. Note the distinction between self-similar and strongly self-similar: the former is a property 
of the overall shape of a single Tile type, the latter is a property of the shapes of individual regions 
in a set of Tile types. 



In the following section we present an active tile assembly system based on a well-known aperiodic 
tiling of the plane obtained by a substitution tiling of an L-shaped tile (FigjHj). The substitution 
process consists of repeatedly "inflating" the L-shaped tile to obtain a larger L-shape. The larger 
L-shape is in fact produced by joining four smaller L-shaped tiles together: every subsequent tiling 
level is obtained from the previous one by substituting each L-shaped tile with four newly assembled 
L-shaped tiles. This is indicated in Figj5] where each level of the substitution is outlined with a 
different color. The tiling of the plane generated in this process is aperiodic^. 



Figure 5. Self-Similar L-Shapc Tiling 

We provide a recursive set of L-shape Tile types and an ATAS which generates them. Then, 
we discuss the recursive properties of the system and prove that it yields an aperiodic tiling of the 
plane. 

3.1. Active Tile Assembly System for the L-Shape Tiling. The premise in the unit tile 
construction that we present here is illustrated in Fig. [6} We begin by joining three unit tiles 
together to obtain the smallest (level 0) L-shape. Then, one of two types of borders is built around 
each L-shape (yellow or red in the figure) and depending on the type of border, the L-shape assumes 
either the center role (yellow) or an outside role (red) in the next level of the assembly process. 
Each subsequent level of the assembly is formed by three outside L-shapes and one center L-shape. 
To proceed to the next level, a border is built once again around each of the assembled L-shapes 
and the process is repeated. We treat the two bordered versions of the smallest L-shape as two 
distinct Tile types in addition to the unbordered type. Fig. [6] outlines three iterations of the process 
yielding levels 0, 1 and 2 of each Tile type. 

We show that the active tile assembly system which produces the desired effect is (Tq, s, 2), with 
To and s given in Fig. [7j Formally, the tiles are given in Tables [5] and [6] in the Appendix. 



3. L-SHAPE 
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Figure 6. Outline of three iterations of the L-shape unit tile construction. A yellow border 
represents the center role and a red one represents the outside role. Each successive level 
unbordered shape (left column) is composed of three red-bordered shapes (center column) 
and one yellow-bordered shape (right column) of the preceding level. 

All of the tiles in the figure (and the two tables) are presented modulo rotation, so, for example, 
the tile G3, given by: 

= ({2},0) 

t_, = ({-66},0) 
t_. = (0,{-4}) 

■A = {Ky} 

^ ~ {^0 ^' 55_y} 
stands to also represent the tiles given by: 



t+y 


= ({2},0) 


t+y 


= ({-66},0) 


t+y 


= (0,{-4}) 


t+x 


= ({-66}, 0) 


t+x 


= (0,{-4}) 


t+x 


= ({i},0) 




= (0,{-4}) 


t-y 


= ({i},0) 


t-y 


= ({2},0) 


t—x 


= ({i},0) 


t—x 


= ({2},0) 


t—x 


= ({-66},0) 


A 




A 




A 


= {4:^} 


§ 


= {^o"^! ^^-y' ^^+x} 


S 


= {2(1^^1 55^^, 55^^} 


S 


= {2o ^1 55^y, 55^^} 



For the remainder of this paper, we adopt the following notational convention: if r is a unit 
tile in Figure [7j then rjy refers to r exactly as shown and te,ts,t\y refer respectively to 1, 2, and 
3 counter-clockwise 90° rotations of Tjy. Thus, for example, GSat is G3 as in the first case listed 
and the three rotations are 03^;, G3s, and G3w respectively. Thus, there are 28, up to rotation, 
equivalence classes encompassing a total of 112 unit tiles required for this tiling (only 28 DNA tiles 
in application, however, since physical DNA tiles are free to rotate in the plane). 
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Active Label 



Inactive Label 
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Initiation Signal 



Activation Signal 



Figure 7. The 28 rotational representatives for the set Tg of unit tiles for the L-shape 

tiling. The strength function s is shown on the left: s(c) = 1 for c = 1, 2, 3, 4, 5 and s(c) = 2 

for c = 11, 22, 33, 44, 55, 66, 77, 88. Note that a black edge (s = 0) represents {A, I) = (0, 0), 

. J.e., the aitsenefl of any labels, ^nd is not jB, label in. itself. , ii, r i j i i n on oo 
As shown mTig. [7J laBeis I, 2, J, 4, ana o are assigned a strength oi 1, and labels 11, 22, 66, 



44, 55, 66, 77, and 88 are assigned a strength of 2 (see Def. |2.l[ ). Note that these labels are to be 
interpreted only as strings of symbols and do not hold any significance as numbers. 

Furthermore, the temperature parameter in (To,s,2) set to "2" means that the sum of the 
strengths of all bonds which would have to be cut in order to separate a tile assembly into two 
parts would have to be at least 2. A consequence of this is that at the initial stage, only unit tiles 
with complementary labels of strength 2 can bind. 

The initial assembly is shown in Fig. l8l under Time 1. The four assemblies on the left are 



complete (see Def. 2.9) assemblies as soon as their unit tiles bind because there is no signaling 
to be transmitted; the assembly X2-X1 requires one iteration of the tile modification function, as 
shown, to become a complete assembly; and the remaining two assemblies become complete after 
twcj^ iterations of /. Thus, Ti is the set To plus the seven (modulo rotation) two-tile assemblies 
added at Time 1. 

Remark. Note that the word "Time" here counts the assembly stages. The assembly figures depict 
only the new tile assemblies added at each stage. The tile assemblies which can interact at a given 
stage i (or "Time" i) are the tile assemblies in Tj and are given by the union of all tile assemblies 
shown in the preceding stages. 

At the next stage, Time 2 (Fig. [s]), the only new bonds that can form are between XI and the 
tile assembly X2-X3 or between X3 and the tile assembly X2-X1. The completed assembly in either 
case is the first L-shape: X1-X2-X3 (in the first case, it takes two iterations of /, in the second only 
one). It can be easily verified that no other additions can occur at this stage: that would require 
side by side labels on one of the two-tile assemblies to be matched with those of another two-tile 
assembly and none of the seven assemblies can be paired in that way. 



Further stage by stage assembly is given in detail in Section 3.3 For now, we state the following 
key properties about the system. 

Theorem 3.1. The active tile assembly system (To,s,2) as given in Fig. yields a recursive set 
Q of self-similar L-shape Tile types with a strongly self-similar subset. 



The first iteration activates labels 2 and 4 for assemblies C0-G3 and C0-G4 respectively and removes ineffectual 
signal patliways. Tfie second iteration removes the ineffectual labels. 
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Figure 8. Time 1 shows the assembhes added to Ti. Time 2 shows the assembhes added 

„j to To, the first L-shape assembhes .in the. cpnstKuction^ • ,i ^ i- ^- rrrrK i 

We give tne construction oi a recursive set for ine system m trie next section (hection 6.2] and 



show that it is, in fact, produced by the system in Section 3.3 



Theorem 3.2. All of the assemblies produced by the L-shape ATAS are components of the L-shape 



Tile types in set from Theorem 3.1 



The above theorem claims that the system produces only structures that take part in the assembly 
of some level of a Tile type in 0. This is proved in Section [3. 3[ 



Theorem 3.3. The L-shape ATAS induces an aperiodic tiling of the plane. 



We prove this in Section 3.4 



3.2. A recursive set of L-shape Tile types. In the following discussion we write tiT2 to repre- 
sent the only possible assembly formed by the two tiles in the ATAS (To, s, 2): e.g., CIatBItv refers 
to the assembly formed by CI and Bl in Figure [8] and X27vX37v to the assembly formed by X2 and 
X3 in the same figure (note that since there is only one way for these unit tiles to form an assembly, 
there is no ambiguity in notation). Also, for any assembly T, T + (n, v) will mean a translation of 
the class representative of T by the vector (u, f). 

We now proceed to define = {Ti, . . . , T12}, a set of 12 (L-shaped) Tile types for (Tq, s, 2), to 
prove Theorem |3.1[ 
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We let the four T, for i = 1 mod 3 represent 
the four possible rotations of the "unbordered" 
Tile types, for z = 2 mod 3 - the "AO-type bor- 
dered" (or edge) Tile types, and for i = mod 3 
- the "EO-type bordered" (or center) Tile types. 
The "AO" and "EO" in the name refer to the 
unit tiles AO and EO in Figure [7] which initiate 
the construction of the corresponding border in 
each case, discussed in Section |3.3[ In the fol- 
lowing sections, we discuss only Ti, T2, T3 in full 
detail as the other Tile types are equivalent up 
to rotation. 

The number of regions for all Tj is specified 

as: 

{4, i = 1 mod 3 
12, i = 2 mod 3 
17, i = mod 3 

Thus, Ti consists of 4 regions, T2 consists of 
12, T3 of 17, and so on. Recall that the jth 
region of Tj, denoted Rij, a configuration, is 
defined as the union of assembly instances r^^ , 
the subregions. We give the number of a re- 
gion's subregions in terms of level £ in Table [T] 
as r]ij{i). 

The type-specific descriptions are grouped 
into three sections: in each, we give Tj(0) ex- 
plicitly for all values of i and also give the re- 
gion and subregion descriptions, accompanied 
by a diagram, for a rotational representative of 





i 


mod 3 


1 


2 
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1 


1 




2 


1 


1 


1 




3 


1 


3-2^-4 


3-2^-4 




4 


1 


1 


1 




5 


- 


3 • 2^+1 - 6 


1 




6 


- 


1 


3 • 2^+1 - 7 




7 


- 


1 


1 




8 


- 


3 • 2^+1 - 6 


1 


3 


9 


- 


1 


1 




10 


- 


1 


1 




11 


- 


3-2^-3 


3 • 2^+1 - 7 




12 


- 


1 


1 




13 






1 




14 






1 




15 






1 




16 






3-2^-3 




17 






1 



Table 1. Table of values for 'qij{£). 



four types. 



ri(o) r4(o) Trio) rio(o) 



1 




X3n 

X27V Xl^v 


1 

y 


X3e 


XIe 
X2e 


1 

y 


XI5 


X25 
X35 


1 

y 


X2vi/ 


X3vi^ 




1 





1 







1 







1 



X X X X 

Table 2. Table oi £ — assemblies for the unbordered Tile types. 



3.2.1. The unbordered Tile types. In this section, we look at Ti,T4,Tt, and Tiq. The initial - level 
- configurations ri(0), T4(0), r7(0), Tio(O) are given in Table [2| Recall that a Tile type at a given 
level is an assembly representative (see end of Section 2.2), a map from the integer lattice into the 
set of active tiles. As shown in Fig.|8j these initial assemblies are obtained in 'J2, the second stage of 
the ATAS (Tq, s, 2). Refer to Figure [9] for a diagram of the four regions of Ti (^) , T4 (^) , (^) , Tio (^) 
for £ > 1. These regions for Ti are given as follows: 
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(a) 



(b) 



Figure 9. (a) The regions of Ti{£),T4{i),Tr{£), and Tio(^) for £ > 1. (b) Level n L-Shape 
Tile Ti{n) + (1,1) in as a complete tile assembly. Tiles marked with represent 
suppressed rows and columns consisting of identical tiles. The only signaling present is on 
the three tiles as indicated by the red and yellow arrows, i.e., signals of the form 1q, 4q 
(see Section 3.3). The strength 1 labels -1,-2,-3,-5 are marked with red, blue, green, and 
light purple edges respectively. The only strength 2 label, -55 at (3 • 2" — 2, 3 • 2"+^ — 4), is 
marked with a dark purple edge. No active or inactive labels other than the ones indicated 
are present. 



^1,1 w 


= rU^) 


= Tu{i- 


- 1 


) + (0,3- 


2^-2) 


^1,2 (^) 


= ri2ii) 


= T2{i- 


1) 


+ (0,0) 




^1,3 (^) 


= rUi) 


= n{£- 


1) 


+ (3 • 2^ 


-2,0) 




= rUi) 


= m- 


1) 


+ (3 • 


-i_ 1,3.2^-1-1 



Recall that Tii,T2,T5 are AO-type bordered shapes and that Ts is EO-type bordered (compare to 
the construction outlined in Figure [6| center and bottom shapes in the left column). Figure [ojshow^s 
the outline of a general Ti(n) Tile type. The Tile types Ti,T-j, and Tio are equivalent to Ti up to 
rotation, so their region descriptions are analogous to the above and are thus omitted. 



3.2.2. AO-type bordered Tile types. We continue with T2,T^,T%, and Tn. The level instances 
T2(0), T5(0), T8(0), Tii(O) are given in Table [3} Refer to Figure 10 for a diagram of the regions of 
T2{t) for i > 1. All regions of T2 with the exception of i?2,i consist of unit tile subregions forming 
a border around the i?2,i region, which is the Tile type Ti (see Figures 
given as follows: 



10 



11). The regions are 
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Table 3. Table oi £ = configurations for tlie AO-bordered Tile types. 
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Figure 10. The regions of T^^i) and outlines of T.^{i),T^{i) and Txi{i) for i>\. The 
regions marked in color consist of a number of subregions that increases with level. Note 
also that the two regions marked in red and blue (i?2,5 and i?2,8 respectively) each consist 
of two nonadjacent sets of subregions. 
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Figure 11. Assembly of T2(n) Tile type: the AO- type border assembles around the Level 
n L-Shape Ti{n) + (1, 1) (i.e., Ti(n) positioned with bottom left corner tile at (1, 1)). 



3-2^+1-6 

R2d^)= U ^5 
k=l 

= rff~^{i) + (0,2) (this shift by 2 units reflects the 
skipping of the single unit tile region iis.e, Fig. flo|, r^t^{() = r|5(^) + (0, 1), /c / 3 • 2^ - 3; 



where rUi) = A2n + (0, 2), rlf-^ 



3-2*+i-6 



R2Ai)= U 4i 



k=l 



where r^8(£) = A3Ar + (1,0), r||'~2(£) = r^|'~3(£) + (2,0) (this shift by 2 units reflects the 
skipping of the single unit tile region i?2,9, Fig. 10), '''2 8^(^) = '^f 8(-^) + (1' 0)' ^ / 3 • 2^ — 3; and 



3-2*-3 



R2M^)= U "^l 
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A:=l 
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where rl -^iii) = Alw + (3 • 2^+^ - 3, 1) and r^^l{£) = r^^iii) + (0, 1). Refer to Figure |lT| for an 
outline of the T2{n) type as a tile assembly. The constructions for T^jTg, and Tn are, once again, 
analogous to the above and are omitted here. 

As we shall see in Section 3.3 both in the case of AO-type and EO-type bordered assemblies, the 
actual construction of the Tile type at a level "starts" with the unbordered assembly in region 1 
binding to a unit tile (AO or EO) in region 2, followed by the region 3 being filled in one tile at a 
time. In our case, the filling of the regions is sequential by design. But from the perspective of the 
recursive Tile type set, it does not matter how or in what order the regions are filled, only that it 
is possible to obtain structures with all regions filled in the indicated way. 
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Table 4. Table of £ configi 
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for the EO-bordered Tile types. 




Figure 12. The regions of T3(£) and outlines of Tii{l) ,Ts){i) and Tx2{t) for i > I. The 
regions marked in color consist of a number of subregions that increases with level. Just 
as in Figure 10 the two regions marked in red and blue {R^fi and -R3.11 respectively) each 
consist of two nonadjacent sets of subregions. 
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3.2.3. EO-type bordered Tile types. We continue with TsjTgjTg, and T12. The construction is es- 
sentiahy identical to that of the Tile types T2,T^,T^, and Tn, but using different unit tiles in the 
border. The initial structures r3(0), T6(0), r9(0), ri2(0) are given in Table [1} Refer to Fi gure 12 
for a diagram of the regions of T-i{() for £ > 1. They are explicitly given as follows: 
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where rl^^i) = EIe + (0,2), ~^(^) = rf^ ~'^{i) + (0,3) (this shift by 3 units reflects the 

skipping of the two single unit tile regions R-sj and -^3,81 Fig. 12), f^'e^i^) = '^f + (Oj 1)' ^ / 
3-2^-4; 

3-2«+i-7 



^3,11 (^) = U 



' 3,11 



k=l 



r. 



3-2^-3 
3,11 



{I) + (3, 0) (this shift by 3 units reflects the 



where rl,^{£) = Els + (1,0), r|;fr'W 
skipping of the two single unit tile regions i?3,i2 and i?3,i3, Fig. 12), = '"t ii(^) + (l) 0), ^ / 

3-2^-3; and 



3-2*-3 



R3M^)= U ^-3, 
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13 



shows the outhne 



where ^gl^) = Elw + (3 • 2^+^ - 3, 1) and r^^jgW = rl i^{^) + (0, 1). Figure 
of the general 73(71) assembly. As with the previous cases, the construction of regions for T^^Tg, 
and T12 follows analogously and is omitted here. 




Figure 13. Assembly of T-i{n) Tile type: the EO-type border assembles around the Level 
n L-Shape Ti{n) + (1, 1) (i.e., Ti{n) positioned with bottom left corner tile at (1, 1)). 

This concludes the construction of the recursive set of Tile types for Theorem |3.1[ What remains 
to be shown for the proof of this theorem is that these structures are, in fact, obtainable from the 
ATAS (To, s, 2). We do this in the next section. 

We finish this section by pointing out that Tj G is self-similar for all i, since the geometrical 
shape of the union of all regions does not change for any i. Furthermore, the set containing any of 
the Tile types Ti, T4, T7, Tiq is a strongly self-similar subset of 6: every region in one of these four 
Tile types consists of a single tile type in 0, but every Tile type in is geometrically similar to 
every other type, so the conclusion follows. 

We point out that the entire set is not strongly self-similar because AO- and EO-type bordered 



Tile types contain regions that stretch, e.g., i?2,3 (see Figure 10), and that aren't geometrically 



similar to any Tile type in (refer to Section 2.5 for the definitions of self-similarity) 



3.3. Self- Assembly. Our goal now is to show that the structures we defined with Tile types do 



in fact assemble in the system (Tq, s, 2), finishing the proof of Theorem 3.1 The construction will 
proceed by induction and show that structures that do not take part in the assembly of Tile types 
in cannot be produced, proving Theorem 3.2 The set of unit tiles is shown in Figure [7| (also 
Tables [5] and [1] in the Appendix) . 
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Figure 14. L-shape assembly, Time 3-6. The completed constructions depicted are the 
only constructions which can be added to for corresponding i. Beginning with this figure, 
only the initial and completed assemblies are shown in all figures, unless indicated otherwise: 
intermediate steps of the tile modification function are suppressed. 

For the base case of the induction we need to obtain ri(l) from the level Tile types. We have 
already seen that XlArX2ArX3Ar assembles in Figure [s] (Time 2), which is Ti{0). We proceed to 
demonstrate the assembly of T2(0) and r3(0) (all other Tj(0) are rotated copies of one of the three 
level types mentioned and assemble in the same way). 

We have already established that no assemblies other than the ones in Fig. [8] could be formed from 
pairing unit tile and double assemblies. Now, observe that all active labels of Ti(0) are strength 1 
except for the label -55 at the top, which has strength 2. All of the strength 1 labels are "negative" 
labels. A "negative" label can only bind to a "positive" label, so two "positive" labels would have 
to be present on the same side of a two-tile assembly in order to bind to Ti(0). The only positive 
labels appropriately positioned on the existing assemblies are the +2 labels on C0-G3, but they 
do not match the labels on Ti(0). Therefore, the only possible binding that can occur is with the 
single strength 2 label, —55. 

We can see that the complementary +55 label appears in both EO and AO unit tiles in Figure [7| 
The two different tiles allow the creation of two different borders around Ti (0) depending on which 
one of EO or AO attaches to Ti (0) . This local non-determinism allows us to give the L-shapes one 
of two different roles (the center and edge roles, respectively) in the next assembly stage. Note that 
border initiation is the only case when two tile choices are available for binding at a single site. 
The assembly process is shown in Fig. 14 Time 3-5 shows the only possible binding events and the 
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Figure 15. L-shape assembly. Time 7 (bottom left), Time 8 (top left). Time 9 (right). 

Again, these are the only new constructions that can be added to the active tile set. 
resulting complete tile assemblies. Note that signals are transferred from AO to CI in Time 4. No 

structures other than the ones shown can be constructed from earlier assemblies because there are 

no assemblies which contain two appropriately positioned strength 1 "positive" labels to match two 

of the "negative" ones, except in the newly formed corners. So, the construction can only proceed 

by matching labels in these corners (which may also contain a "negative" and a "positive" strength 

1 label). An activated strength 2 label does not occur until the formation of the next level Tile 
type. 

Throughout the rest of the L-shape assembly, we use the fact that the unbordered Tile types 
only have sequences of strength 1 "negative" labels along their sides except for the single strength 

2 label -55 which initiates the creation of the AO or EO border. So, only this label of strength 2 or 
the new corners formed by two negative (or, sometimes, negative and positive) labels can serve as 
attachment points for other tiles. 

Note that in Figure 14, Time 6, both the AO-type assembly and the EO-type assembly may be 
augmented by two different assemblies: either the unit tiles B2 and G3 respectively, or the two-tile 
assemblies formed at Time 1, B2 and C2 assembly or G3 and CO assembly. However, in the former 
case. Time 7 would produce the latter case. Thus, the only new assemblies in Time 7 (Figure 15) 
are the ones from Time 6 assemblies that were created with the two-tile assemblies. 

This process, when continued, yields the desired r2(0) and Ts{0) assemblies (shown in Tables [S] 
and|4]) at Time 9 in Figure 15 The rotated Tile types (not shown) are obtained simultaneously in 
a similar manner. 

We can now assemble ri(l). In the T^iO) assembly (and in the corresponding rotations) we see 
two perpendicular label 4's (yellow) (negative on the EO tile and positive on the Gl tile), which 



24 



N. JONOSKAt.i, D. KARPENKO*'t>i 



can be matched by the Tii(O) assembly (a rotation of T2(0)), as it has a +4 on the AO and a —4 
on the D3 tile. When these attach (see Fig, 16, Time 10), a signal from Gl activates the label — 1 
(red) on C2 and C3, and a (yellow) signal from D3 is sent through Gl to G3, activating the —4 
label in the completed Time 10 (not shown). Now the T^iO) portion of the Tii(0)T3(0) assembly 
in Time 10 has active -4 and +4 labels in tiles G3 and G2 respectively, allowing T2(0) to attach at 
Time 11 (Fig. 16). A signal from G2 activates label —2 (blue) on C2 and C3, and label —5 on CI 
in the completed Time 11 (not shown). The —4 label (yellow) on G4 is activated similarly to G3. 
Finally, in Fig. 17 (left), the third AO Tile type, 75(0), can attach to the -4 label in G4 and and 
+4 label in F3. A signal from F3 activates label —3 (green) on C2 and C3, and label —5 on CI. 



Also, a signal travels from G4 to AO, then back to G4 and eventually to CI on the Tii(O) (Fig. 17 
top left), activating the —55 label (Fig. 17, right). In this way, we obtain Ti(l) in Figure 17 



Note that starting with the Time 10 stage, we also obtain incomplete L-shape constructions from 
completed AO-type assemblies binding to the EO-type assemblies with only a partially assembled 
border (for example, to the EO-type shape like in Time 7, Fig 15). These are not depicted as they 
ultimately lead to constructions that result in the desired L-shapes. An example of the partial 
border mixed assemblies can be seen in Figures 32|p6 in the Appendix and also in |12j . 



CO + - EO 



I . I 



T^IO (start) 



1 ,1 



- .±. 



AO |+_Q3 



T=1 1 (start) 



1 .1 



Figure 16. Time 10-11, assembling the level 2 shape. The AO-type assembly (rii(O)) 
binds via AO and D3 to EO and Gl of 13(0). The red signal from Gl is transmitted to 
C2 and C3, activating the corresponding labels. A (yellow) signal from D3 is transmitted 
through Gl and activates the -4 label in G3 which allows an AO- type assembly (T2(0)) to 
bind at G3 and G2 sites. Broken transmission signals and disconnected inactive labels are 
cleaned up. In the figure on right, a process similar to the previous takes place, resulting in 
the activation of designated labels on the CI, C2, and C3 tiles in the AO-type and the G4 
tile in the EO-type. 

We can now proceed with the inductive step. Suppose we begin with a shape like in Figure[9j that 
is, with Ti{n) for some n > 1 (clearly, Ti(l) satisfies the diagram) and suppose we have T2(n — 1) 
and T^{n — 1) (and all corresponding rotations). We first show that from Ti(n) we can obtain both 
T2(n) and T3(n). This is demonstrated in Figures 11 and 13 Note that Ti{n) (Figure [9]) contains 
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only negative strength 1 labels around its borders except for the strength 2 label, -55, in the top 
right corner. So, Ti{n) cannot interact with any assemblies except for the unit tiles AO and EO, 



yielding two possible border initiations, shown in Figures 11 and [13} Upon close examination, it is 
evident that beginning with the AO or EO tile in the corner, it is possible to construct the borders 
one (or two, where possible) unit tile at a time exactly as indicated. 

Moreover, observe that the exposed edges of the growing assembly in either figure never contain 
"+" strength 1 labels except for the +4 (yellow) labels on AO (Fig. 11), and Fl, F2, and F3 unit 
tiles (Fig. 13 Note that although each of these three tiles appears twice, the +4 label is activated 
only the second instance of each, after receiving a signal from the C4 tiles). In Figure 13 , an exposed 
-4 complementary label first appears on EO, parallel to the L-shape outline. Additional active -4 
labels on this assembly are exposed only one at a time and are perpendicular to the L-shape border 
while building the EO border. Here we note that it is geometrically impossible for any assembly 
to bind simultaneously to a perpendicular positive or negative 4 label and a parallel positive or 
negative 4 label on another assembly. The reason is as follows: if the labels appear on the same side 
of the L-shape border, e.g., parallel +4 on Fl (or -4 on EO) and perpendicular -4 on El (during the 
border building process, just as one of these unit tile attaches), the matching assembly would have 
to have a parallel -4 (or -|-4) and a perpendicular +4 exposed, but no exposed perpendicular -|-4 
labels ever appear on any assembly. If the labels appear on different sides, the matching assembly 
would have to round two or more corners, whereas the existing L-shape assemblies are only capable 
of rounding one corner. So, only pairs of parallel labels 4 in Figure 13 remain to be considered. 
Because none of the labels 4 appear on the same side, the EO-bordered assembly is geometrically 
prevented from binding to its own type as it assembles. On the other hand, Figure [TT] shows that 
the AO-bordered figure has a perpendicular -|-4 label on AO and acquires a perpendicular -4 label 
on D3 at the last stage of its assembly, allowing it to bind to the +4 and -4 labels on adjacent 
sides of the EO-bordered assembly, specifically, to the EO and Fl (the only correctly positioned 
complements) . 




Figure 17. Time 12 completes the assembly of ri(l) unbordered L-shape from four level 
bordered shapes. 

What this shows is that no level of either assembly (EO or AO bordered) can interact with another 
level except for an AO Tile type to bind to an EO Tile type of the same level: first at EO and Fl, then 
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at E2 and F2, and finally at E2 and F3 (this sequence of binding events is enforced by signaling - E2 
is only activated when the binding at EO and Fl takes place). Observe that the next border cannot 
be initiated until the unbordered shape is completely assembled, producing Ti(n + 1). The reader 
is invited to compare the tile assemblies in Figures [TT] and [13] to the region outlines in Figures [T0| 



and 12 respectively. (72(1) and T3(l) are shown as an example in Figure 35 in the Appendix). 

The position coordinates for the tiles in the L-shapes are derived from the fact that Ti(0), the 
assembly obtained at Time 2 (see Figure [S]), measures 1x2 tiles on the short and the long sides 
respectively, and if Ti{l) measures L x 2L tiles, then the bordered versions T2{1) and T^^i) each 
measure (L + 1) x (2L + 2) tiles. 

Let Si denote the length of the short side of Ti{£). Since Ti{i + 1) is obtained by combining four 
of the bordered shapes of the previous level as in Figure^), we get the following relation: 



si = 2se-i + 2, for £ > 1 
= 3-2^-2, for £ > 

where the second equality can be easily verified using the first: so = 1 = 3 • 2^^ — 2, so if = 3-2^ — 2, 
^ > 0, it follows that se+i = 2(3 • 2^ - 2) + 2 = 3 • 2^+^ -4 + 2 = 3- 2^+^ - 2 as required. 

Based on the above, we can also compute the precise stage at which a Tile type at a given 
level assembles. For example Ti(0) assembles at stage 2 and for £ > 1, Ti{i) assembles at stage 
9 • 2^+1 -6i- 18. 

So far, we have shown that the Tile types of the previous section do result from the self-assembly 
process of this tiling system. We have also shown why this tiling system does not produce any 



assemblies that are not components of the Tile types in 0. This concludes the proof of Theorems 3.1 
andIO 



As a final note here, we include a side by side comparison of the assembled Ti(l) and Ti(2) 



structures in Figure 18 The assembly process from Ti(l) up to Ti(2) (stage 42) is shown explicitly 



in the Appendix and can also be found in [T2j. 

3.4. Aperiodicity. In this section we give the proof of Theorem |3.3[ First, we show that the 
L-shape Tile types provide a tiling of the plane by considering the Ti (unbordered) Tile type. We 
note that for each i, the corresponding L-shape can be divided into three square regions, and that 
each of the square regions in an £ + 1 L-shape covers more than 4 times the area covered by a square 
region of the level i L-shape. Thus, since these squares grow without bound, given any sized square 
portion of the plane, there exists some m such that the level m L-shape covers it. As m — )• oo, the 
L-shape assembly instances provide a tiling of the planej^ 



For the proof of aperiodicity, refer to Figure 18 and Figures [TT] and 13 Notice that for a given level 



£ L-shape Ti(£), there is a unique outermost C2 marked with a —2 label and a unique outermost C3 
marked with a —2 label. These two unit tiles, furthermore, are separated by a region of A3's, where 
the exact number of the A3's indicates the level i. This uniqueness is guaranteed by construction: 
every level i L-shape is formed by an EO-type bordered level i — 1 L-shape, T3{£ — 1) in the center 
that is adjacent to three AO-type bordered level £ — 1 L-shapes, T2{£ — 1),T5(£ — l),Tii(£ — 1) 
(Fig. [9]). Each of the AO- type shapes contains a single C2 and a single C3 as part of its border. The 
label counter-clockwise adjacent to the +22 (dark blue) label on C2 and to the +33 (dark green) 



label on C3 (see Fig, 18 ) is determined by the AO-type level £ — 1 shape's position relative to the 



level £ — 1 EO-type shape. The three possible orientations yield three possible label values (-1,-2, 



^This is, in fact, a consequence of Konig's lemma. 
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Figure 18. Ti (1) and Ti (2) L-Shapes side by side. Note the positioning of corner markings 
(red labels -1, blue labels -2, green labels -3, and the purple -55 label) and observe that the 
rest of the outside tile edges are filled with the -5 label. 

or -3). Consequently, the particular configuration of C2 and C3 carrying the —2 label is indeed 
unique within the level i L-shape. 

Fix an L-shape tiling of the entire plane, ^. We consider Ti{£) and Tt{£). Observe that T'i'{i) 
(refer to Fig. [9^)) appears inside Ti{£ + 2) (it is nested inside T8{£) that is inside Tii{£ + 1)) by 
construction. The unique placement of C2 and C3 followed by the fixed number of A3's in ri(^) 
and Tili) guarantees that no translation of an L-shape tiling of the plane is possible by a vector 
with either component less than 3-2^ — 2 (the length of the short side of a level i L-shape Ti 
or Tj) for any i (Ti and Tj together cover the four possible orientations of the C2 C3 segment). 
Otherwise, supposing this is not the case for some £ >2, consider the plane tiling induced by the 
L-shape assembly instances. Then a C2 C3 border configuration as described above for the level 
i L-shape would have to appear more than once inside that level £ L-shape Ti (if at least one of 
the vector coordinates is positive) or T7 (if both vector coordinates are negative) - and this, we 
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have shown, cannot happen. Since this is true for every level it follows that the tiling ^ is not 
invariant under any translation. 



4. Conclusion 

In this paper, we presented an abstract self-assembly model for active DNA-based tiles - tiles with 
capacity to receive and transmit binding site activation signals. We also gave a definition of recursive 
assembly and self-similarity based on this abstract model, which, to the best of our knowledge, is 
the first attempt to describe recursion in self-assembly. The aperiodic L-shape assembly system 
that we described illustrates not only these new concepts of active tiles and recursive tilings but 
also the power of signaling in the self-assembly process. It was the signaling capability of the tiles 
that ensured that no extraneous structures could arise and enforced the order of assembly without 
the necessity of staging. 

We point out that the method used for the creation of the L-shape Tiles can be used to create 



a variety of aperiodic patterns: for example the square pattern in Figure 19 The principle is 
to use local nondeterminism to create the desired hierarchy. At each level, each instance of the 
shape is assigned a role in the next level by obtaining a partial border, selected at random, which 
grows in length with proportion to the size of the shape. This growth in length is what guarantees 
non-periodicity as increasingly longer border lines consisting of one kind of tile types will be found 
further away from a given tile, disallowing translational symmetry. Essentially, the entire identity 
of the larger tile is contained in its outside border, so the borders mediate all binding and signaling 
which occurs among the shapes. 




Figure 19. A non-periodic square tiling utilizing the border technique; here, four different 
border types - red, yellow, green, and blue - determine a particular square's position in the 
larger shape. 
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As mentioned in the introduction, this model can be adjusted to accommodate different as- 
sumptions to better approximate the physical self-assembly process. For example, the alternate 
description of an active tile model made in [1^ includes the ability to deactivate what we here 
have called "edge labels" and allowing tiles to break away from a constructed assembly. In fact, 
adding a deactivation capacity to the model we have presented is straightforward. If desired, the 
triple (r, yi,S) could be expanded to include a set of deactivation signals and a set of deactivation 
pathways, becoming a 5-tuple. The tile modification function would have to be expanded to re- 
move deactivation signals when they are used and to remove labels that have been deactivated in 
a manner completely analogous to activation signal processing. However, after any deactivation, 
it would be necessary to check whether the assembly is still ^-stable - if not, the simulation of 
the remaining signal transmission would need to be continued on separate assemblies (which, if 
not stable themselves, would need to be broken up as well). The definition of active tile sets (the 
stage hierarchy of sets) would also have to be modified to include assemblies which are obtained in 
this manner. We did not include deactivation in the model presented here in order to simplify the 
exposition in this first introduction to the model and to recursion in self-assembly. 

Likewise, one could account for the variability in signal transmission rates, as the authors in 
|19) do, by allowing tiles to "wait" before transmitting a signal. Our assembly model is rather 
simplified in the sense that we assume that all signal transmissions occur simultaneously and that 
all signaling across a single assembly is completed before it can bind to another one. In solution, all 
these processes of binding and signal transmission are occurring in parallel and at variable rates; 
thus, a more accurate model of physical reality would allow an assembly at any point in its signal 
transmission cycle to bind with other assemblies, which means for any assembly a we would have 
to add assemblies /*(«) for every z > (where / is the tile modification function) to the active 
supertile set of the relevant stage. An even more accurate model could account for a different rates 
of signal transmission by writing instead of Cg whenever a tile receives an activation signal from 
another tile; this could represent a pending signal. A tile containing such a signal could replace the 
X with a at any iteration of the tile modification function, in which case it would act the same 
way as the already defined activation signals. 

We note that neither of the possible extensions that we mentioned would affect the assembly 
process in the L-shape ATAS that we described in this paper since the design of the tiles does not 
allow two Tile types at a given level to join before both of them are completed (with respect to 
the tile modification function, i.e., signal transmission). Thus, in our case, the simpler model is 
sufficient for an accurate presentation of the given recursion. 

References 

[1] G. Aggarwal, M.H. Goldwasser, M.Y. Kao, and R.T. Schweller. Complexities for Generalized Models of Self- 
Assembly. Proceedings of the Fifteenth annual ACM-SIAM symposium on Discrete algorithms, p. 889, 2004. 

[2] R.D. Barish, P.W.K. Rothemund, and E. Winfree. Two Computational Primitives for Algorithmic Self-Assembly: 
Copying and Counting. NanoLetters, 5:2586-2592, 2005. 

[3] B. Chakraborty, N. Jonoska, N.C. Seeman. A Programmable Transducer Self-Assembled from DNA. Chem. Sci., 
3:168-176, 2012. 

[4] E.D. Demaine, M.L. Demaine, S.P. Fekete, M. Ishaque, E. Rafalin, R.T. Schweller, and D.L. Souvaine. Staged 
Self-Assembly: Nanomanufacture of Arbitrary Shapes with 0(1) Glues. Natural Computing, 7(3), 2008. 

[5] E. Dolzhenko, N. Jonoska, and N.C. Seeman. Transducer Generated Arrays of Robotic Nano-Arms. Natural 
Computing, 9:437-455, 2010. 

[6] T.J. Fu and N.C. Seeman. DNA Double Crossover Structures. Biochemistry, 32:3211-3220, 1993. 



30 



N. JONOSKAt.i, D. KARPENKO*'t'i 



[7] K. Fujibayashi, R. Hariadi, S.-H. Park, E. Winfree, and S. Murata. Toward Reliable Algorithmic Self- Assembly 

of DNA Tiles: a Fixed- Width Cellular Automaton Pattern. NanoLetters, 8(7):1791-1797, 2008. 
[8] C. Goodman-Strauss. Matching Rules and Substitution Tilings. The Annals of Mathematics, 147:181223, 1998. 
[9] B. Griinbaum and G.C. Shephard. Tilings and Patterns. Freeman, 1987. 
[10] H. Gu, J. Chao, S.-J. Xiao, and N.C. Seeman. A Proximity-Based Programmable DNA Nanoscale Assembly 

Line. Nature, 465(7295) :202-205, 2010. 
[11] N.R. Kallenbach, R. I. Ma, and N.C. Seeman. An Immobile Nucleic Acid Junction Constructed from Oligonu- 
cleotides Nature, 305:829-831, 1983. 
[12] D. Karpenko. Self-assembly of Self-similar Structures by Active Tiles. Master's Thesis. University of South 
Florida. 2012. 

[13] T.H. LaBean, H. Yan, J. Kopatsch, F. Liu, E. Winfree, J.H. Reif, and N.C. Seeman. The Construction, Analysis, 
Ligation and Self-Assembly of DNA Triple Crossover Complexes. J. Am. Chem. Soc., 122:1848-1860, 2000. 

[14] K. Lund, A.J. Manzo, N. Dabby, N. Michelotti, A.Johnson-Buck, J. Nangreave, S. Taylor, R. Pel, M.N. Sto- 
janovic, N.G. Walter, E. Winfree, and H. Yan. Molecular Robots Guided by Prescriptive Landscapes. Nature, 
465:206-210, 2010. 

[15] U. Majumder, T.H. LaBean, and J.H. Reif. Activatable Tiles for Compact, Robust Programmable Assembly and 

other Applications. LNCS 4848:15-25, 2008. 
[16] C. Mao, T.H. LaBean, J.H. Reif, and N.C. Seeman. Logical Computation using Algorithmic Self-Assembly of 

DNA Triple-Crossover Molecules. Nature, 407:493-496, 2000. 
[17] T. Omabegho, R. Sha, and N.C. Seeman. A Bipedal DNA Brownian Motor with Coordinated Legs. Science, 

324(5923):67, 2009. 

[18] J. Padilla, W. Liu, N.C. Seeman. Hierarchical Self Assembly of Patterns from the Robinson tilings: DNA Tile 
Design in an Enhanced Tile Assembly Model. Natural Computing, online first, DOI: 10.1007/sll047-011-9268-7, 
2011. 

[19] J. PadiUa, M.J. Patitz, R. Pena, R.T. Schweller, N.C. Seeman, R. Sheline, S.M. Summers, and X. Zhong. 
Asynchronous Signal Passing for Tile Self- Assembly: Fuel Efficient Computation and Efficient Assembly of 
Shapes. 

Available on Arxiv: |http: //eirxiv. org/pdf /1202 . 5012vl .pdf] 

[20] L. Qian and E. Winfree. A Simple DNA Gate Motif for Synthesizing Large-Scale Circuits. DNA Computing, 
LNCS, 5347:70-89, 2009. 

[21] L. Qian and E. Winfree. Scaling Up Digital Circuit Computation with DNA Strand Displacement Cascades. 

Science, 332(6034):1196-1201, 2011. 
[22] R.M. Robinson. Undecidability and Nonperiodicity for Tilings of the Plane. Inventiones Mathematicae, 12(3):77- 

209, 1971. 

[23] P.W.K. Rothemund, N. Papadakis, and E. Winfree. Algorithmic Self-Assembly of DNA Sierpinski Triangles. 

PLoS Biology,2{12):e424, 2004. 
[24] P.W.K. Rothemund. Folding DNA to Create Nanoscale Shapes and Patterns. Nature, 440(7082) :297-302, 2006. 
[25] G. Seelig, D. Soloveichik, D. Y Zhang, and E. Winfree. Enzyme-Free Nucleic Acid Logic Circuits. Science, 

314(5805):1585, 2006. 

[26] N.C. Seeman. DNA Nicks and Nodes and Nanotechnology. NanoLetters, 1:22-26, 2001. 

[27] W.B. Sherman and N.C. Seeman. A Precisely Controlled DNA Bipedal Walking Device. NanoLetters, 4:1203- 
1207, 2004. 

[28] J.-S. Shin and N.A. Pierce. A Synthetic DNA Walker for Molecular Transport. J. Am. Chem. Soc, 126:10834- 
10835, 2004. 

[29] H. Wang. Notes on a Class of Tiling Problems. Fundam. Math., 82:295-305, 1975. 

[30] Y. Wang, J.E. Mueller, B. Kemper, and N.C. Seeman. The Assembly and Characterization of 5-arm and 6-arm 

DNA Junctions. Biochemistry, 30:5667-5674, 1991. 
[31] E. Winfree. Algorithmic Self- Assembly of DNA. Ph.D. Thesis. California Institute of Technology. 1998. 
[32] E. Winfree, F. Liu, L. A. Wenzler and N.C. Seeman. Design and Self- Assembly of Two-Dimensional DNA Crystals. 

Nature, 394:539-544, 1998. 

[33] H. Yan, X. Zhang, Z. Shen, and N.C. Seeman. Directed Nucleation Assembly of DNA Tile Complexes for 
Barcode-Patterned Lattices. Proc. Nat. Acad. Sci. 100:8103-8108, 2003. 



ACTIVE TILE SELF-ASSEMBLY, SELF-SIMILAR STRUCTURES AND RECURSION 



31 



[34] P. Yin, H. Yan, X.G. Daniell, A.J. Turberfield, and J.H. Reif. A Unidirectional DNA Walker that Moves Au- 
tonomously Along a Track. Angewandte Chemie International Edition, 43:4906-4911, 2004. 

[35] B. Yurke, A.J. Turberfield, A.P. Mills, F.C. Simmel, and J.L. Neumann. A DNA-fuelled Molecular Machine 
Made of DNA. Nature, 406(6796) :605-608, 2000. 

[36] H. Zhong W. Liu, R. Wang, and N.C Seeman. Crystalline Two-Dimensional DNA Origami Arrays. Angew. 
Chemie, 50:264-267, 2011. 



32 



N. JONOSKAt>i, D. KARPENKO*'t.i 

Appendix 



Tile Index: Tile Sides with Active and Inactive Label Sets 



Tile 


t+y 




t^y 


t—x 


XI 


(0,0) 


({-3},0) 


({-3},0) 


({77}, 0) 
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({-88}, 0) 
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/ r o1 rA\ 
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/Tool rh\ 
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G\ 
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/ r r 1 rA \ 
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({-i},0) 


/ r Ai rfi\ 
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G2 


/ r r 1 (K\ 
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({-2},0) 


G3 
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(0,{-4}) 


G4 


({3},0) 


({-66}, 0) 


(0,{-4}) 
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Fl 


(0,{4}) 
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/ r 1 1 rA\ 

({i},0) 


(0,{-5, -44}) 


F2 
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(0,{-5,-44}) 


(0, {4}) 


F3 


({3},0) 


(0,{-44}) 


(0,{4}) 


({4},0) 


EO 


({-4},0) 
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/fit r ^ r<^ \ 

(0,{-4, -66}) 


El 
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({4} , 0) 
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E2 
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//"A rA\ 

(0,0) 


CI 


(0,0) 
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63 
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Al 


(0,0) 


({i},0) 


({5},0) 


({-i},0) 


A2 


({2},0) 


({5},0) 


({-2},0) 


({-5},0) 


AS 


({5},0) 


({-3},0) 


({-5},0) 


({3},0) 



Table 5. List of unit tile rotation class representatives for the L-Shape Tiling: tile sides 
and labels. 
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Tile Index: Transmission and Activation Signal Sets 


Tile 
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X3 
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AO 


{-'-0 ^ -1 ^0 ^' ^^+!e} 




Al 
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A2 
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A3 
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Table 6. List of unit tile rotation class representatives for the L-Shape Tiling: signaling. 
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Figure 20. Time 13-14, initiating the AO and EO border formation. The only two possible 
binding events that can occur is the level 2 L-shape binding to an AO or an EO. Any other 
events are disallowed either by geometry or by complementarity and strength: observe that 
all of the labels on the outside of the L-shape are "negative" strength 1. Two "positive" 
strength 1 labels would have to match between it and another structure in order for a bond 
to form. The reader may verify that none of the structures which could be produced up to 
this point have two positive labels in the correct configuration. 
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Figure 21. Time 15-16, continuing AO-EO border formation. The assembly of the borders 
continues. Since Bl may already be attached to CI, both the single and the double addition 
structures are produced in Time 16; but the former may produce only the latter at the next 
step. Note also that the positive-positive corners required to fit into the negative-negative 
adjacencies are only ever present on single and double tiles, never on any larger structures 
(most tiles have only two positive labels and they use these to attach to growing assemblies, 
thus no positive labels remain for binding; the ones that do not follow the pattern are G2, 
G3, and Fl, F2, and F3, but these contain positive 4 labels and there is never more than 
one such label on a fixed side of the assembly, so their geometry forces them to bind only 
with AO- type assemblies of the same level, as seen in Fig 32|[37 1 . 
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Figure 22. Time 17-18, the assembly continues. Note that in 17, the gold label "-44" on 
Fl was activated by a signal in the previous time step, which allowed a C4 tile to attach in 
this step. 
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Figure 23. Time 19-20. 
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Figure 24. When Fl attaches, it receives the signals 4 (yeUow) and 5 (hght purple) 
originally sent by C4 and transmitted by the El's. 
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Figure 25. Time 23-24. A similar situation arises with B2-C2 as did with Bl-Cl. Two 
assembhes are produced: one containing only B2 and one containing the B2-C2 pair; the 
former produces the latter at the next time step; so both evolve into the same shapes, but 
one does so faster. 
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Figure 26. Time 25-26. The signal originally sent by E2 is received by F2 and activates 
label -44 allowing C4 to attach. 
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Figure 27. Time 27-28. 
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Figure 28. Time 29-30. Notice that the second F2 receives the same signals as the second 
F2, activating its +4 and +5 labels and allowing E2 to attach at the next step. 
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Figure 29. Time 31-32. Just like Bl-Cl and B2-C2, C3 can either be already attached 
to B3, or become attached immediately after. 
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Figure 30. Time 33-34. 
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Figure 31. Time 35-36. The AO- type border is completed at Time 36. 
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Figure 32. Time 37. Two things can occur. Either an AO-type assembly binds to a 
partially complete EO-type assembly as in figure on the left, or with even less of the border 
complete (not shown), or an AO- type assembly may gain another border tile (right). 
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Figure 33. Time 38. If the AO- type assembly binds, a signal from D3 is transmitted to 
the nearest E2, activating its -4 label and allowing it to bind to another AO-type bind if the 
second F2 on the EO-type border is present (left). The existing assemblies may also simply 
gain another border tile (right). 




Figure 34. Note that the final AO-type assembly cannot attach until the border on the 
EO-typc assembly is completed, so at this step only border additions may be gained. 
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Figure 35. T2(l) and r3(l) (left and right, respectively). 
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Figure 36. Time 40 (top two assemblies) and 41 (bottom assembly). We see that in 
41 only one new kind of assembly can be produced, because the variants of the partially 
assembled L-shapes converge to the same shape as more binding events occur. 
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Figure 37. Time 42. The first level 3 shape is assembled. Note the similarities between 
this and the level 2 shape. Both are given side by side in Fig 18 



